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Abstract 

We evaluate the pion-nucleon and the pion-Delta sigma terms by employing the method of 
quantum chromodynamics (QCD) sum rules. The obtained value of the pion-nucleon sigma term 
is compatible with the larger values already anticipated by the recent calculations. It is also found 
that the pion-Delta sigma term is as large as the pion-nucleon sigma term. 
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I. INTRODUCTION 



The meson-baryon sigma terms are important for hadron physics as they provide a mea- 
sure of chiral-symmetry breaking and the scalar quark condensate inside the baryon. In 
particular, the pion-nucleon sigma term has received much attention and has been exten- 
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llj and various other approaches 
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221 ]. It is equivalent to the value of the scalar form-factor 
m(N{p')\üu + dd\N{p)) = aN{k)v{p')v{p), 



m 



m u + m d 



(1) 



at zero momentum transfer viz. a^ = cr(0), where v{p, s) is the Dirac spinor for the nucleon, 
m is the average quark mass and k = {p' — p) 2 . The pion-nucleon sigma term is also defined 
via the Feynman-Hellmann theorem as 

dm¡sf 



CTjv 



q=u,d 



m 



dm n 



m{N\üu + dd\N), 



(2) 



where m^ is the nucleon mass. In the chiral limit, the sigma term vanishes; however, the 
fact that pions and the u- and the c¿-quarks are not actually massless implies that a^ ^ 0. 
The pion-nucleon sigma term cannot be measured directly and there is no consensus on its 



20| 



value. Various estimates in the literature range from 18 ± 5 MeV [3| to 74 ± 12 MeV 

As for the pion-nucleon sigma term, the pion-Delta sigma term is of recent interest. The 
Delta resonance, A(1232), plays an important role in pion-nucleon scattering away from 



24 



25 



261 ] . It is close in mass to the 



threshold and in loop calculations of %PT 0, l2i 
nucleon and it couples strongly to nucleons, photons and pions. In this respect, the value 
of the pion-Delta sigma term is crucial in relation to Delta-nucleon mass splitting as the 
sigma term provides a measure of the shift in the hadron mass away from the chiral limit. 
It is equivalent to the value of the Delta scalar form-factor at zero momentum transfer viz. 
°a = o"a(0), where the latter is defined as 



m(A(p', s')\üu + dd\A{p, s)) = -v,{p', s')\g^a A {k) + p w p»F T {k)\v v {p, s). 



(3) 



Here, a A {k) and F T {k) are the scalar and tensor form-factors, respectively, and v^{p,s) is 
the Rarita-Schwinger spin-vector of the Delta, with the spin projection s. The minus sign 
on the right-hand side (RHS) is conventional like in the case of the free Delta Lagrangian. 



Analogously to Eq. ([2]), the pion-Delta sigma term is defined via the Feynman-Hellmann 
theorem as 

mí¡^ = -m(A(s)\üu + dd\A(s')), (4) 



0"A 



q=u,o 



and it corresponds to the change in Delta mass with the quark mass. Recently, some 
model-dependent approaches have been used to obtain <ta : chiral-quark model produces 
o"a = 32 ±3 MeV 17J], which is in agreement with the value from the so-called configuration 
space quark model, oa = 30 ±2 MeV 27J]. A preliminary chiral perturbation theory analysis 
of the lattice data gives o"a = 20.6 MeV 28| . A solution of the Faddeev equation produces 
a relatively larger value as oa — 50 MeV [22]. 

Our aim in this work is to calculate the pion-nucleon and the pion-Delta sigma terms 
using the external-field QCD sum rules (QCDSR), which are a powerful tool to extract 
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32]. In this 



qualitative and quantitative information about hadron properties 
framework, one starts with a correlation function that is constructed in terms of hadron 
interpolating fields. On the theoretical side, the correlation function is calculated using the 
Operator Product Expansion (OPE) in the Euclidian region. This correlation function is 
matched with an Ansatz that is introduced in terms of hadronic degrees of freedom on the 
phenomenological side. The matching provides a determination of hadronic parameters like 
baryon masses, magnetic moments, coupling constants of hadrons, and so on. To determine 
the value of the sigma terms, we evaluate the vacuum-to-vacuum transition matrix element 
of two nucleon and Delta interpolating fields in an external isoscalar-scalar field. Our anal- 
ysis here follows closely the one in Ref. [l5[, where was calculated using this method and 
<jn = 36 ± 5 MeV was obtained. We depart from this prior analysis by 1) incorporating 
the dimension-7 operators 2) improving the QCDSR analysis with up-to-date values of the 
vacuum parameters, in particular the vacuum susceptibilities and 3) treating the uncertain- 
ties more systematically by employing a Monte-Carlo analysis. Moreover, we report on the 
value of the pion-Delta sigma term in QCDSR by following the same approach. 

Our paper is organized as follows: In Section [TTl we present the formulation of QCDSR 
and construct the relevant sum rules. We give the numerical analysis of the sum rules and 
discuss the results in Section IIIII Finally, we arrive at our conclusions in Section [IVl 
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II. THE DERIVATION OF THE SUM RULES 



In the external-field QCDSR method, one starts with the correlation function of the 
baryon interpolating fields in the presence of an external constant isoscalar-scalar field S, 
defined by 



i j d 4 x é p - x (0 \T[ VN (x)rj N {0)}\0) s = U(p) + SU s (p) + 0(S 2 ), 

i j d'x é** (o W A (xW A (o)}\ o) 5 = n""(p) + sn s »( P ) + o(s 2 ), 

where r¡ N and rj A are the nucleon and the Delta interpolating fields, respectively: 

VN = tabc [u^C-f^U b ] ^l^d c , 
Va = e abc [u^C'J^Ub] U c . 



(5) 



(6) 



Here a, b, c are the color indices, T denotes transposition and C = Í7 2 7°. For the interpolat- 
ing field of nucleon, there are two independent local operators, but the one in Eq. (¡6]) is the 



optimum choice for the lowest-lying positive parity nucleon (see e.g. Ref 33] for a discussion 
on negative-parity baryons in QCDSR). Tl(p) and U^^p) are the correlation functions when 
the external field is absent and correspond to the functions that are used to determine the 
baryon masses. The second terms in Eq. (j5J) represent the linear response of the correlators 
to a small external scalar-field S, which is computed with an additional term to the QCD 
Lagrangian, which is 

AC = -S g s q [ü(x) u(x) + d(x) d(x)] . (7) 

Here, S represents the external scalar-field and gjj is associated with the coupling of the 
external scalar-field to the quark. The external scalar-field contributes to the correlation 
functions in Eq. (¡3¡) in two ways: first, it directly couples to the quark field in the baryon 
currents and second, it modifies the condensates by polarizing the QCD vacuum. In the 
presence of an external scalar-field there are no correlators that break the Lorentz invariance, 
like (qa^q) which appears in the case of an external electromagnetic-field F^ u . However, 
the correlators already existing in the vacuum are modified by the external field, viz. 

(qq)s = (qq) -xS(qq), 

(8) 

(g c qa ■ Gq) s = (g c qv ■ Gq) - XcS(g c q(T ■ Gq), 



4 



T[r) N (x)r¡ N (0)} 



2ie abc e a'b'c'r 



2ie abc e a'b'c' 



where x an d Xg are the susceptibilities corresponding to the quark and the quark-gluon 
mixed condensates, respectively, and the coupling of the external scalar-field to the quark 
is simply taken as g q = 1. 
At the quark level, we have 

'Tr{S b u b ' {x)^[S?(x)\ T C^YtrfSf(xW>K, 

Tr{SÍ b '(x)YC[Sf(x)} T C^}Sf(x) (9) 

+ 2S b u b ' (x)YC[ST\x)] T CYS c :\x) 

To calculate the Wilson coefficients, we need the quark propagator S q in the presence of 
the external scalar-field, which is given in Ref. 15j. Using this quark propagator, one can 
compute the correlation function ITs(g) and U s u (q), which can be brought into the form 

Tl s (p) = U N (p 2 )^ + W N (p 2 ), 

(p) = u A ( P 2 ) <ri> + W A (p 2 ) gr + ■ ■ ■ , 

where the ellipsis represents the Lorentz-Dirac structures other than g^ v and g^vj). 

Note that one can obtain the sum rules at different Lorentz structures. Here, we choose 
to work with the sum rules at the structures p 1 and g^ v p 1 for the nucleon and the Delta, 
respectively, where the latter is completely contributed by the Delta baryons with J = § 



(see e.g. Ref. 
as 



32J for details). The OPE sides of the sum rules linear in S are then obtained 



TÍ N (P 2 ) 
TÍa( P 2 ) 
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(11) 



3p 2 2p 2 9p 4 

In the above equations, we have defined the quark condensate a = —(2n) 2 (qq), and the 
quark-gluon-mixed condensate (qg c a ■ Gq) = ml(qq) with the QCD coupling-constant 
squared g 2 c = 4iia s . The four-quark condensate is parameterized as ((qq) 2 ) = ^(qq) 2 . We 
note that the last term of IlAr(p 2 ), which is associated with the dimension-7 operators has 



been neglected in Ref. [15]. 



The analyticity of the correlation function allows us to write the phenomenological side 
of the sum rules in terms of a double-dispersion relation of the form 



ReU B (p 2 



1 



OO POO 



lmTL B (p) 



(si -p 2 )(s 2 -p 2 [ 



dsi ds 2 . 



(12) 



The ground-state hadron contribution is singled out by utilizing the zero-width approxima- 
tion, where the hadronic contributions from the Breit-Wigner form to the imaginary part of 
the correlation function is proportional to the 5-function: 

lmU N (p) = 7r 2 c^i -m^)^ -m^KOlr^ 

+ vr 2 5(si - m%)5(s 2 - m N *)(0\r) N \N(p))(N(p)\S(üu + dd)\N*(p))(N*(p)\rj N \0). 

Here, the second term is associated with the transitions to higher nucleon states. We then 
expresses the correlation function for the nucleon as a sharp resonance plus a continuum 
after Borel transformation: 

n K (M^(A> K f + c K ^) L^l + l£^ e -^ it , (14 ) 

by using the definition in Eq.(¡2]). The correlation function for the Delta is similarly expressed 
as 

In Eqs. (Tbil) and (ÍT5|) . the matrix elements of the currents r¡ N and 77 A between the vacuum 
and the hadron states are defined as 

(0\r¡ N \N(p,s)) = X N v(p,s), 

(16) 

(0|í&|A(p,s)) = A A iAp,s), 

respectively, for the nucleon and the Delta, where X N and Aa are the residues. We also make 
use of the Rarita-Schwinger spin-sum, which is 

£ v»(p, s)v»(p, s) = - (gT ~ \l»l v ~ P \ - flÍl ~ ^ff) + ™a). (17) 

\ 3 3 m A 3 m A / 

We have included the single-pole contributions with the factors C N and C¿\, which corre- 
spond to the transitions to higher baryon states. These transition terms are not properly 
suppressed after the Borel transformation and should be included on the phenomenological 
side. 

The QCD sum rules are obtained by matching the OPE sides with the hadronic sides 
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and applying the Borel transformation. The resulting sum rules are 

-aM*E - ÍM! x «a 2 L*/» - ^ M 2 a L~^ 27 + ^ a 2 ( X + Xc ) L^ 

6 (18) 



( XlmN^- + C N M 2 + 5w 2 L ^ e <rt,-*V*A 

\ m 2 / 



e 



3a M 4 E L 16/27 - ^! x/í a 2 L 28 / 27 - ^M 2 aL 2 / 27 + ^ a 2 ( X + Xg) L 14 ' 27 

3 2 2 9 (19) 

= ( ~\lm A ^ +C A M 2 + —w* 5w 2 L ^¿m\-^)IM-\ ¿^Jtf 
\ m 5 / 

where M is the Borel mass and we have defined A| = 32tt 4 A| with B = N, A. The 
continuum contributions are included via the factor 

E n = l-(l + x + ... + ^)e~ x , (20) 

with x = w 2 /M 2 , where u> is the continuum threshold. In the sum rules above, the third 
terms on the RHS of the sum rules denote the contributions that come from the response 
of the continuum threshold to the external field. Here, 5w 2 represents the variation of the 
continuum threshold and the coefficient is calculated bydifferentiating the continuum parts 



3J | with respect to the quark mass. 



of the chiral-even nucleon and Delta mass sum rules 

These terms are suppressed by the factor e -( w ~- m B)/ M2 as compared to the single-pole terms, 



however, should be included on the phenomenological side if 5w 2 is large (see Ref. 35[] for a 
detailed explanation of this term). The corrections that come from the anomalous dimen- 
sions of various operators are included with the factors L = log(M 2 / Aq CD ) / \og(fi 2 / Aq CD ) , 
where ¡x = 500 MeV is the renormalization scale and Aq C d ís the QCD scale parameter. 

III. ANALYSIS OF THE SUM RULES 

We determine the uncertainties in the extracted parameters via the Monte-Carlo based 
analysis introduced in Ref. 



361 ] . In this analysis, randomly selected, Gaussianly distributed 
sets are generated from the uncertainties in the QCD input parameters. Here we use a = 
0.52 ± 0.05 GeV 3 , b = (g 2 c G 2 ) = 1.2 ± 0.6 GeV 4 , m 2 = 0.72 ± 0.08 GeV 2 , and A QCD = 
0.15 ± 0.04 GeV. The factorization violation in the four-quark operator is searched via the 
parameter n, where we take k = 2 + 1 and 1 < k < 4; here ((qq) 2 ) > (qq) 2 is assumed via 
the cut-off at 1 (for a discussion on QCD parameters see e.g. Ref. (36J). The value of the 



susceptibility X has been calculated in Ref. 



37J]asx = — 10 ± 1 GeV and has been used 



to calculate the scalar-meson-baryon coupling constants [38]. The susceptibility xg is less 
certain, therefore we take its value equal to the one of Xi however with a larger uncertainty 
Xg — —10 ±3 GeV^ 1 . We use 10 4 such configurations from which the uncertainty estimates 
in the extracted parameters are obtained using a fit of the LHS of the sum rules to the RHS. 
For m, we make use of the Gell-Mann-Oakes-Renner relation which is 

2m(qq) = -m 2 Jl (21) 

where m n = 138 MeV is the pion mass and f n = 93 MeV is the pion decay constant. 
We use the following chiral-odd nucleon and Delta mass sum rules for normalization of the 



sigma-term sum rules [39j, |40|] : 



aM 4 -^ab= | m)V e -< /M2 , (22) 

-aE^L^M* - -E m 2 aL 2/27 M 2 - -abL w/27 = S mAe -™Í/« 2 , (23) 
3 3 18 2 



which have been found to be more reliable than the chiral-even sum rules |36|, |41[. The 
Monte-Carlo analysis of the sum rules fíT8|) and fíT9l) is performed by first fitting the mass 
sum rules to obtain the pole residues Ag and these residue values are used in the sigma-term 
sum rules for each corresponding parameter set. We find that the nucleon mass sum rule 
in Eq. (1221) fails to resolve the pole from the continuum. In order to proceed and reduce 
the uncertainties in the final results as much as possible, we fix the nucleon mass at its 
experimental value as mjv = 0.94 GeV and the continuum threshold at w = 1.5 GeV. This 
produces A^ = 1.60 ± 0.18 GeV 6 . It is also found that the mass sum rule (|23|) somewhat 



overestimates the value of the Delta mass 41|. Therefore, instead of fixing the mass at 



its experimental value, we make a two parameter fit including m¿\ and Aa by fixing the 
continuum threshold at w = 1.7 GeV, which gives mA = 1-45 ± 0.05 GeV and A A = 
5.20 ±0.66 GeV 6 . 

The valid Borel regions are determined so that the highest-dimensional operator (HDO) 
contributes no more than about 10% to the OPE side which gives the lower limit on the 
valid Borel region. The upper limit is determined using a criterion such that the continuum 
contributions are less than 50% of the phenomenological side. It is relevant to point out 
that the dominant contributions to the OPE sides of the sum rules in (Í18l) and ( ÍT9l) come 
from the terms that involve the susceptibilities i.e. the second and the fourth terms on 
the LHS. This leads to a suppression of the continuum contributions, which are involved 
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FIG. 1: The continuum-plus-transition contributions as well as those of continuum change for the 
sum rules (ÜSJ) and (ÍTÜj) . 

by the first terms on the LHS only. In Figure [U we plot the continuum-plus-transition 
contributions as well as those of continuum change for the sum rule (11 8p . In the Borel 
window 0.9 GeV< M < 1.3 GeV, which we determine as the fiducial region, the relative 
double-pole-contribution is higher than 50%, while the continuum contributes less than 
about 10%. We also observe that the effects of the continuum variation with the external 
field are minor. For the sum rule (Íl9l . we are not able to find a Borel window according 
to the above criterion: the double-pole contribution is less than 50% in the region where 
the HDO starts to contribute less than 10%. As can be observed in Figure [U in the Borel 
window 1.3 GeV< M < 1.5 GeV which we take as the fiducial region for the sum rule 
analysis, contribution of the continuum-plus-transition to the RHS of the sum rule is about 
60%, while the HDO contributes less than 10% to the OPE side. 

To demonstrate how well the sum rules and the fitting work, we first arrange the sum 
rules in the subtracted form 

m (24) 
m 

where represents the OPE-minus-excited and minus-continuum-change contributions. 
As the RHS appears as a straight line with this form, the linearity of the LHS gives an 
indication of the OPE convergence and the quality of the continuum model. In Figure [2, 
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we present the fit of the sum rules ( ÍTBI) and (ÍTÜI) . respectively, using the average values of 
the QCD and the obtained fit parameters. The error bars at the two ends correspond to the 
uncertainties in the QCD parameters. 

Since the continuum contributions in the sum rules ífTSl) and ÍÍTD"]) are suppressed as 
compared to the total phenomenological side, it becomes difficult to extract information 
about the continuum threshold from the fit. Therefore, we have assumed that the continuum 
thresholds are equivalent to those for the mass sum rules. The variation of the continuum 
threshold Sw 2 can also be obtained from the fit. However, instead of taking this as a free 
parameter, we proceed with a generous assumption that the continuum threshold changes 
by 25% with the external field viz. 5w 2 = w/4. Then, a two parameter fit of the sum rule 
(ÍTgj) including the a N and C N from a consideration of 10 4 parameter sets produces 

a N = 53 ± 24 MeV, (25) 
and a fit of the the sum rule ( ÍTÜi) including the <ta and Ca gives 

a A = 54 ± 25 MeV. (26) 
Eq. f[2"5]) suggests a 60% enhancement of the pion-nucleon sigma term as compared to the 



value obtained in Ref. 15] , which is 36 ±5 MeV. The large errors of about 50% in our results 
(as compared to typical 30% in QCDSR) mainly stem from the uncertainties in the residue 
values and the vacuum susceptibilities. We also find that the dimension-7 terms contribute 
by 20% to the OPE side. As emphasized above, the continuum effects are suppressed, 
therefore a change in the value of the continuum threshold e.g. by 10%, leads to a negligible 
change in the final values. One can also obtain a ratio of the pion-nucleon sigma term to 
the pion-Delta sigma term by dividing the corresponding values of the two sigma terms at 
each QCD parameter set. The analysis of the final distribution obtained from this method 
produces 

o-jv/íta = 1-06 ±0.16, (27) 

which is in accordance with the observation that the pion-Delta sigma term is as large as 
the pion-nucleon sigma term. 
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FIG. 2: The subtracted form of the sum rules (I18p and ([19]). The solid line is the double-pole 
contribution and the dashed-line is the OPE-minus-excited states and minus-continuum-change 
contributions, where we use the average values of the QCD and the obtained fit parameters. The 
error bars at the two ends represent the uncertainties in the QCD parameters. 

IV. CONCLUSIONS AND DISCUSSION 

We have calculated the pion-nucleon and the pion-Delta sigma terms using the external- 
field QCDSR method. Our analysis includes several improvements over the prior work 
in Ref. which enhance the central value of the pion-nucleon sigma term by 60%. We 
observe that inclusion of the dimension-7 operators in the OPE analysis brings a contribution 
of 20% to the final results. A Monte-Carlo analysis of the QCD input parameters provides 
a more systematic treatment of the errors. Although the errors in the final results are large 
(which could be best improved by a more precise determination of the residues), a number of 
qualitative features are evident. The value we have obtained for the pion-nucleon sigma term 
is consistent with the larger values already anticipated by the recent calculations. Moreover, 
our analysis favors a pion-Delta sigma term as large as the pion-nucleon sigma term. 



In Ref. 42J, inconsistencies were pointed out in the leading nonanalytic behavior in m q 
between the treatment of the nucleon mass with the QCDSR and the ^PT description. The 
inconsistency arises because the quark condensate has a chiral behavior of the form 

— m_ In — 

327r 2J2 M Q 



(qq) = (qq) ( 1 - m£ ln ^ + . . . ) , (28) 
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whereas the nucleon mass is known to have a leading nonanalytic term proportional to m' q 2 
in the form 

M N = M N + Aml ^Ar- ml + ..., (29) 

but no m q \nm q term. (Here g¿ is the nucleon axial-vector coupling constant, M is a mass 
parameter of the order of ~ 1 GeV, A is an unknown constant 43| and all quantities with a 



circle over denote the first term in the chiral expansion of that quantity.) These spurious and 
missing nonanalytic terms lead to an intrinsic uncertainty of the order of ~ 100 MeV in the 
QCDSR estimates of M N and a N . It was shown in Ref. ¡44^ with a more general argument 
that the spurious chiral-log contributions to M N originate from virtual pions, which cancel 
if the continuum contribution to the correlator is treated carefully by including the 7r-N 
continuum via soft-pion theorem. Unlike the case of M N , it is not straightforward to prove 
the same for the current sum rule for a N , by explicitly including the 7r-N continuum. The sum 
rule in Eq. (fT8|) takes account of the change of w with m q in a complicated way. As stated 



in Ref. 



421 ]. this m q dependence of w may eliminate the discrepancy between the QCDSR 



and the xPT descriptions. One should, however, bear in mind that the usual continuum 
model is inconsistent with the treatment of the continuum by including the virtual pions 44¡ , 
and difficulties arise with such an arbitrary solution. On the other hand, our approach is 
equivalent to evaluating the derivative in the Feynman-Hellmann theorem using a QCDSR 
estimate of M N , which yields a QCDSR estimate of a N . Such a relation between the sum rule 
of M N and that of a N implies that spurious chiral-log contributions to a N should also cancel 
with the inclusion of the 7T-N continuum. However, we note that, although the spurious 
terms and the resulting uncertainties are removed by a proper treatment of the continuum, 
the problem of the missing m q ^ 2 term, which originates from the chiral expansion of M N 
in Eq. (129p . persists. We know of no satisfactory method for recovering this missing term, 
which is left as an open problem in QCDSR. The lack of this term leads to uncertainties of 
~ 15 MeV and ~ 20 MeV for M N and a N in QCDSR [42]. 

The sigma term is the measure of the contribution of explicit chiral-symmetry breaking in 
the baryon masses. The QCD Hamiltonian consists of the chiral-invariant terms containing 
the gauge couplings of gluons and the chiral-non-invariant quark-mass term. Suppose that 
the chiral non-invariant term is weak and therefore treated perturbatively. Then the sigma 
term is nothing but the contribution of the quark-mass term to the baryon mass. As we find 
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that the pion-Delta sigma term is of similar size with the pion-nucleon sigma term, the quark- 
mass term is concluded to give little contribution to the Delta-nucleon mass difference. This 
in turn means that their mass difference is entirely due to the gluon-quark gauge-coupling 
term. Its main role is to induce the spontaneous chiral-symmetry breaking, giving the quark 
constituent mass, but the perturbative part is known to yield spin-spin (color-magnetic) 
interaction between the quarks. The latter is well-known to cause the Delta-nucleon mass 
difference in the quark model. 

Finally, we note that the central values of the sigma terms obtained in Eqs. (|25|) and (|26|) 
correspond to the 6% and 4% of the physical nucleon- and the Delta-masses, respectively, 
which would indicate the change in these baryon masses if the chiral symmetry is restored. 
We plan to extend the QCDSR analysis to hyperons and calculate the hyperon sigma terms, 
as well, to address the explicit chiral-symmetry breaking in the complete baryon octet 45]. 
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